Gravito-magnetic instabilities in anisotropically expanding fluids 
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Gravitational instabilities in a magnetized Friedman - Robertson - Walker (FRW) Universe, in 
which the magnetic field was assumed to be too weak to destroy the isotropy of the model, are known 
and have been studied in the past. Accordingly, it became evident that the external magnetic field 
disfavors the perturbations' growth, suppressing the corresponding rate by an amount proportional 
to its strength. However, the spatial isotropy of the FRW Universe is not compatible with the 
presence of large-scale magnetic fields. Therefore, in this article we use the general-relativistic (GR) 
version of the (linearized) perturbed magnetohydrodynamic equations with and without resistivity, 
to discuss a generalized Jeans criterion and the potential formation of density condensations within 
a class of homogeneous and anisotropically expanding, self-gravitating, magnetized fluids in curved 
space-time. We find that, for a wide variety of anisotropic cosmological models, gravito-magnetic 
instabilities can lead to sub-horizonal, magnetized condensations. In the non-resistive case, the power 
spectrum of the unstable cosmological perturbations suggests that most of the power is concentrated 
on large scales (small k), very close to the horizon. On the other hand, in a resistive medium, the 
critical wave-numbers so obtained, exhibit a delicate dependence on resistivity, resulting in the 
reduction of the corresponding Jeans lengths to smaller scales (well bellow the horizon) than the 
non-resistive ones, while increasing the range of cosmological models which admit such an instability. 

PACS numbers: 04.25.Nx, 04.40.Nr, 04.20.Jb 



I. INTRODUCTION 

It is known that, the formation of large-scale structures 
(such as galaxies, super-clusters of galaxies and/or super- 
clouds of stars within galaxies) both in cosmological and 
in astrophysical scale is closely related to the concept of 
instability [1], [2]. 

The first attempt towards a theory of galaxy forma- 
tion carried out by Sir James Jeans in the early nine- 
teenths [3], [4]. He assumed that the Universe is filled 
with a self-gravitating non-relativistic fluid, character- 
ized by its mass-density (p), the corresponding pressure 
(p) and a velocity field (u). As a consequence, the re- 
sulting dynamical system is governed by the equation 
of continuity, together with Euler's and Newton's equa- 
tions for the matter-content and the gravitational field, 
respectively [5]. However, Jeans' theory does not apply 
to the formation of galaxies in an expanding Universe, 
because it involves a static medium. Nevertheless, he 
arrived at the conclusion that, density fluctuations with 
wavelengths A greater than a critical value Aj will grow, 
so that the system, eventually, becomes unstable (grav- 
itational instability). As a consequence, any isothermal 
gaseous sphere with a length-scale greater than Aj is 
gravitationally unstable and is going to contract [6]. 

The first satisfactory theory of gravitational instabil- 
ities in an expanding Universe was proposed by Lif- 
shitz [7]. He showed that, in a homogeneous and isotropic 
Universe filled with a perfect fluid, cosmological pertur- 
bations with wave-numbers (fc) below the Jeans' one (fcj) 
grow, not exponentially, but, as a power of the corre- 
sponding scale factor, R(t). Among the various scenar- 
ios of cosmic fluctuations there is also the case where 



gravitation and cosmic expansion are essentially irrele- 
vant. In this context, Weinberg's work [8] is of particular 
interest: For the first time, the notion of an imperfect 
fluid was used and the role of dissipation in the survival 
of protogalaxies was taken into account, along the lines 
of Eckart's formalism [9], which differs from the corre- 
sponding of Landau & Lifshitz [10]. Accordingly, Wein- 
berg found that the protogalactic fluctuations, which be- 
have as ordinary sound waves during the period where 
their matter-content is much lower than the correspond- 
ing Jeans mass (M << Mj), damp in the acoustic phase 
due to the viscosity of the imperfect fluid. 

The study of gravitational instabilities in the presence 
of magnetic fields in an expanding Universe arose as a 
natural generalization [1 1] . Accordingly, solutions for the 
evolution of the mass-density and the magnetic field (H l ) 
perturbations (in the linear regime) have been quested in 
FRW cosmological models, in which, the magnetic field 
was assumed to be too weak to destroy the isotropy of 
the model [12] - [16]. 

Although it has been argued that a uniform magnetic 
field slows down the growth-rate of the unstable modes, 
several authors have studied the influence of cosmological 
magnetic fields (either homogeneous or inhomogeneous) 
to the formation of structures within a cosmic medium 
(e.g. see [17] - [25], for an extensive, though incom- 
plete list). In this context, it was shown [2G] that several 
magnetohydrodynamic processes like pinch-effects, hose- 
instabilities, sausage- and kink-instabilities, may also 
play a significant role to structure formation. These pro- 
cesses correspond to individual harmonics in the eigen- 
function solutions of the wave equation for plasma in- 
stabilities and although they are different, their time- 
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development is similar. 

In the present article, we intend to discuss gravita- 
tional instabilities in a homogeneous and anisotropic cos- 
mological model, in which the large-scale anisotropy is 
due to the presence of a homogeneous magnetic field 
along the ^-direction. 

Our choice is justified by the fact that, mathemati- 
cally speaking, the spatial isotropy of the FRW Universe 
is not compatible with the presence of large-scale mag- 
netic fields. In fact, an anisotropic cosmological model 
can and should be imposed for the treatment of mag- 
netic fields whose coherent length is comparable to the 
horizon length [27]. Therefore, although current observa- 
tions give a strong motivation for the adoption of a FRW 
model, the effects one may lose by neglecting the large- 
scale anisotropy induced by the background magnetic 
field, should be investigated. Not to mention that the 
anisotropy of the unperturbed model facilitates a closer 
study of the coupling between magnetism and geometry. 

The Paper is organized as follows: In Section II, we 
summarize the gcneral-relativistic version of magnetohy- 
drodynamics (MHD) and present its application to the 
Thome's class of homogeneous and anisotropic cosmo- 
logical models [28]. They correspond to a class of semi- 
realistic solutions to the Einstein-Maxwell equations, ad- 
equate to describe a model which departs from isotropy 
along the radiation epoch [29]. In Section III, we perturb 
the corresponding MHD equations to study the evolution 
of linear fluctuations propagating along the i-direction, 
i.e. perpendicularly to the magnetic field. An immediate 
result is the derivation of the corresponding dispersion 
relation, which may reveal the critical values of the per- 
turbations' wave-numbers in order to discuss on their po- 
tential growth (or decay). We find that, for a wide class 
of anisotropic cosmological models, gravito -magnetic in- 
stabilities can lead to sub-horizonal, magnetized conden- 
sations. In this case, the power spectrum of the unsta- 
ble cosmological perturbations suggests that most of the 
power is concentrated on large scales (small fc), close to 
the horizon. These results are in contrast to what has 
been found in the isotropic case, where the magnetic field 
disfavors large-scale matter aggregation, by reducing the 
growth-rate of the corresponding density perturbations 
[11], [13], [14]. In Section IV, we consider the relativis- 
tic MHD equations within a resistive medium. Now, the 
critical wave-numbers so obtained exhibit a delicate de- 
pendence on resistivity, resulting in the reduction of the 
corresponding Jeans lengths to smaller scales than the 
non-resistive ones, while increasing the range of cosmo- 
logical models which admit such an instability. 



II. MAGNETOHYDRO DYNAMICS IN A 
MODEL UNIVERSE 

In the system of geometrical units, where both the ve- 
locity of light (c) and Newton's constant (G) are equal to 



unity, the GR version of the MHD-cquations [18] reads 

(p u a uP = h^(p + ^-), aP + 2^(H 2 6) 

-(H a HP) iaf3 + 2 X (— + <T 2 -u J 2 -u a u a ) (1) 
+f (P + 3p + H 2 ) + 2u a {H a H f3 ), p + {H 2 ). a u a 

supplemented by the equation of motion T a S = 

±u a + X u a + X 0u a + ( P + ^- ) .,pg aP - (H a H ) . p = (2) 
and Maxwell's equations in curved space-time 

h" = « +un- m v + -4- p,» H " ufl ( 3 ) 

6 p + p 

In Eqs (1) - (3), the dot denotes differentiation with re- 
spect to t, Greek indices refer to the four-dimensional 
space-time (in accordance, Latin indices refer to the 
three-dimensional spatial section) and 5% is the Kro- 
neckcr symbol. Furthermore, p is the mass-energy den- 
sity as measured in the rest-frame of the fluid, p is the 
corresponding pressure and are the components of 
the magnetic field. Accordingly, x stands for the com- 
bination x — P + P + H 2 , u a is the fluid's four- velocity 
and h a/3 = g a @ + u a uP is the projection tensor. Finally, 
a^u and lo^ v stand for the shear- and the rotation- ten- 
sor respectively, while 6 = u a a is the expansion param- 
eter, with the semicolon denoting covariant derivative. 
Accordingly, u a = u%vfi . The above system is supple- 
mented by the Raychaudhuri equation [30] 

u% = e+^ + 2(a 2 -u 2 ) + ±(p + 3p + H 2 ) (4) 

Eqs (1) - (4) govern the evolution of magnetized plasma- 
fluids within the context of GR. In what follows, we 
use them to describe the propagation of small-amplitude 
waves (cosmological perturbations) in a homogeneous 
and anisotropic cosmological model, the anisotropy of 
which is due to the presence of an ambient magnetic field 
(frozen into the matter-content) along the ^-direction, 
H = H{t)z. 

The corresponding line-element is written in the 
form [28], [29], [31] 

ds 2 = -dt 2 + A 2 (t){dx 2 + dy 2 ) + W 2 {t)dz 2 (5) 

In Eq (5), the scale factors A(t) and W(t) are solutions to 
the Einstein-Maxwell equations for a magnetized perfect- 
fluid obeying the equation of state 

P = 1P (6) 

Therefore, the various models of this class are param- 
eterized by 7 taking values in the interval 7 G [| , 1). 



3 



Accordingly, A(t) ~ t 1 / 2 and W(t) ~ where I 
(l-7)/(l + 7 ) [28], while 

3-7 



and 



P(7,*) 



and 



2t 2 (l + 7 ) 2 

(l- 7 )i/2(3 7 - i)V2 
2*(l + 7) 



(7) 



(8) 
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FIG. 1: The energy-density of the magnetic field (H /8n) in 
cgs units, at a fixed time (t ~ 10 -21 sec), as a function of 7 
taking values in the range (| , 1). It becomes maximum in a 
cosmological model with 7 = 0.6. 

Notice, that for 7 = | the strength of the magnetic 
field vanishes and the model reduces to the radiation- 
dominated FRW model. As t — > 0, these models evolve 
towards a singularity of infinite density and magnetic 
field. However, for every i < 7 < 1, the ratio 
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(1-7) (3 7 -l) 
16tt (3 - 7) 



(9) 



remains (a non-zero) constant, suggesting that the influ- 
ence of the magnetic field to the evolution of the curved 
space-time is comparable to that of the matter. 

We assume (e.g. see [29]) that the metric (5) can de- 
scribe the early Universe prior to the recombination epoch 
(tree ~ 10 13 sec). At t = t r ee (he. far away from the ini- 
tial singularity) the magnetic field is too weak to have 
any effect on the dynamics of the model, but preserves 
the inherent anisotropy. Therefore, the spatial section of 
(5) always remain axially-symmetric. In fact, it is the 
magnetic field that accelerates the transverse expansion 
(over the xy-plane), while decelerating the longitudinal 
one along the z-direction [28]. 

For an observer comoving with the fluid in the space- 
time (5), we obtain 



(37 -I) 2 
6(l + 7 ) 2 i 2 



(11) 



To study the conditions under which the background (5) 
admits gravitational instabilities or not, we turn to the 
theory of cosmological perturbations. 



III. JEANS LENGTHS IN A MAGNETIZED 
PERFECT FLUID 

For every dynamical system, much can be learnt by 
investigating the possible modes of small-amplitude os- 
cillations or waves. A plasma is physically much more 
complicated than an ideal gas, especially when there is 
an externally applied magnetic field. As a result, a vari- 
ety of small-scale perturbations may appear. 

Accordingly, we perturb Eqs (I) - (4), introducing (to 
the physical variables of the fluid) small-amplitude fluc- 
tuations along the x-axis and neglect all terms of or- 
ders higher or equal than the second. Within the lim- 
its of linear analysis, excitation of cosmological pertur- 
bations in a homogeneous and anisotropic cosmological 
model is basically a kinematic effect, in the sense that 
the self-gravitation of the fluctuations is unimportant 
(e.g. see [6], pp. 501 - 506 and references therein). 
In this case, perturbations' growth arises mostly due 
to their motion in the anisotropic background. There- 
fore, as far as the enhancement of MHD perturbations 
in the anisotropic space-time (5) is concerned, we may 
neglect the first-order corrections of the metric, admit- 
ting the so-called Cowling approximation [32]. In other 
words, in what follows, we treat the MHD perturbations 
as very slowly-varying (low) frequency waves propagat- 
ing in an anisotropically expanding medium, without in- 
teracting with the curved space-time, unless the linear 
regime breaks down. 

To develop the theory of small-amplitude waves in the 
curved space-time (5), we follow the so-called adiabatic 
approximation [1], [34], [35] considering that all the per- 
turbation quantities arc proportional to the exponential 



exp 



ilkx- / n(t') dt' 



(12) 



where, k is the comoving wave-number. In this con- 
text, the (slowly- varying) time-dependent frequency of 
the wave [n(t)] is defined by the eikonal 



n 



n(t') dt' 



(13) 



through the relation n = dQ/dt. Accordingly, we assume 
a wave-like expansion of the perturbation quantities, in 
the form 



Sp = A p e 



t(kx — J t ndt') 



SH 



A H e 



t(kx — J t ndt') 



5u 

Sp ■ 



A, 



a i(kx — /* ndt') 



c 2 s S P 



(14) 
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where c 2 — Sp/Sp = 7 is the speed of sound (in units 
of c). In the anisotropic model under consideration, the 
MHD equations (1) - (4) linearized over the perturbation 
quantities, result in 

Su [2 X rj x - g 33 r§ 3 H 2 + P : Q + in X ] = ik [ 7 <Jp + H SH] 

(15) 

and 



(HSH) [(2 -7) 9 + in] 



1 + 7 



[6(l + ~f)-m}6p (16) 



where are the Christoffel symbols of (5), g^ v are 
the contravariant components of the corresponding met- 
ric tensor, u 2 A = is the (dimensionless) Alfven velocity 

[elf Eq (9)] and we have set p* = p + 

Before attempting to discuss any temporal evolution 
of the perturbation quantities, it is important to trace 
what kind of waveforms are admitted by this system. 
To do so, we have to derive their dispersion relation, 
D(k,n) = [35]. Provided that certain kinds of modes 
(such as acoustic, magnetosonic etc) do exist, they can 
be excited through their interaction with the anisotropic 
space-time. An additional excitation, due to the non- 
zero resistivity, is also possible [26]. We have to point 
out that, although the background quantities depend on 
time, in the search for a dispersion relation, we treat the 
perturbations' amplitudes (Ajs) as constants (at least ini- 
tially). In this way, our search for potential waveforms, 
is not disturbed by the inherent non-linearity introduced 
by the curved background. Nevertheless, once the po- 
tential waveforms are determined, their interaction with 
the curved space-time in the presence of an external mag- 
netic field, implies that, for t > 0, the time-dependence of 
their amplitudes is a priori expected [elf Eq (34), below]. 

The combination of Eqs (15) and (16) results in the 
dispersion relation 



D r (k,n) +iDi(k,n) = 
the real part of which is written in the form 
22 - 9 7 ' 



D r (k, n) = — n 
+ k 2 

where, we have set 
16 



9(2-7) 



Jx — R3 — 



2-7 



(17) 



(18) 



Ri 



Ji = -(l + 7) + \-(l + 3ci)(l + ci)p (19) 



and 



i?i = -p(33 7 + 29) 9 2 + 4<t 2 + 2p (1 + 2 7 ) + ^-H 2 , 
R3 = if [(18 7 - 122)# 2 + 15(1 + 7 M (20) 



The vanishing of Eq (18) leads to a critical wave-number, 
for which the cosmological perturbations under consider- 
ation are stable (oscillating); namely, n 2 > 



k 2 > fc 2 7 (7,*) 
5u = 



- Ji + R 3 + 



■Ri 



(21) 



As a consequence, for k 2 < fc 2 7 , the small-scale fluctua- 
tions of the cosmic-medium become unstable, acquiring 
an imaginary frequency [n 2 (t) < 0]. In other words, they 
are no longer oscillating, but evolve exponentially with 
time (Jeans-like instability). 

In the limiting case where 7 = 5 (FRW radiation 
model), we have H = = u 2 A , a 2 = 0, 9 7^ and Eq 
(18) takes the simpler form obtained in [16]. In this case, 
a Jeans instability arises at wave-numbers k 2 < k 2 ^^, t), 
where 



k 2 rj {\t)c 2 s 



(l + C 2 )(l + 3 C 2 )p+y(l + C 2 )fl 2 (22) 



On the other hand, for 7 — > 1 (Zel'dovich ultra-stiff mat- 
ter) we obtain H = = u\, while a 2 ^ 0, 9 ^ and Eq 
(18) also yields the corresponding result of [16]; namely, 



k 2 j (l,t)c 2 s = (l + C 2 J(l + 3c 2 Jp+U^9 2 + 2o- 2 



(23) 



Finally, in the extreme case of a static model, i.e. 9 = 
and a = 0, we recover the well-known Jackson's solu- 
tion [36]. Accordingly, we conclude that, the critical 
wave- number (21) possesses the correct asymptotic be- 
havior in every limiting case. 

The rhs of Eq (21) is a function of time, as well as of 
7. In particular, 



fis(7)/Qa(7) 



(7 + 1) 2 * 2 
where, for all 7 in the range [^ , 1), 

Q 3 (7) = 36 7 (3 -7) (2 -7) 



(24) 



+ 18 (I-7) (37-I) >0 



(25) 



and 



ffe (7) = 2(1 + 7 )(3 - 7 )(2- 7 )(277 2 + 36 7 + 161) 

+ (1 - 7) (37 - 1)(2 - 7)(-15 7 2 + 174 7 - 931) 
+ 3(1 - 7) (37 - 1)(9 7 2 - IO67 - 99) > (26) 

i.e. they are (positive) polynomials of 7, of the third- 
and the fifth-order, respectively (see Fig. 2). We observe 
that, for every i < 7 < 1, we have 



1 



and therefore 



K.j(l,t) > A, 7 (-,i) 



(27) 



(28) 
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In other words, the more we depart from a FRW ra- 
diation model, the larger is the scale at which gravito- 
magnetic instabilities can become prominent. Departure 
from 7=5) suggests that the initial fortification of the 
magnetic-field strength - for 1 < 7 < 0.6 - is compen- 
sated by a subsequent decrease - for 0.6 < 7 < 1 (in 
connection sec Fig. 1). However, in any case, the corre- 
sponding length-scale grows - at least up to 7 < 0.85, af- 
ter which it remains practically constant (Fig. 2). There- 
fore, we may conclude that, any non-vanishing ambient 
magnetic field disfavors small-scale (sub-horizonal) con- 
densations within a non-resistive plasma-fluid. This re- 
sult is compatible to what has already been found in the 
isotropic case (e.g. see [13], [14]). 



in geometrical units or 




FIG. 2: The critical wave- number A; 2 , (in cgs units) at a fixed 
time (t ~ 10 -21 sec), as a function of 7 taking values in the 
range [| , 1). 

Since k rj ~ the entirety of critical wave-numbers 
in the form (24) vanish only at t — ► 00. Hence, in 
the curved space-time (5) we may always find a length- 
scale involving a gravito '-magnetic instability. However, 
the MHD perturbations considered, correspond to low- 
frequency waves and, hence, these results should be valid 
(mainly) at the late stages of the radiation epoch (prob- 
ably around recombination). For every t within the ra- 
diation era, we may use Eqs (24) - (26) to estimate the 
corresponding critical length-scales (see also Fig. 2). In 
accordance, 



(2.5 



4.7) x - 

' t 



(29) 



(in geometrical units) where, the higher values corre- 
spond to 7 in the vicinity of ^ (in particular, 7 < 0.5), 
while, the lower ones to the rest of the range. The cor- 
responding Jeans length reads 



A rj = - (1-3 - 2.5) x t 

k r j 



(30) 



A r 



(4 - 7.5) x 10 10 x t 



(31) 



in cgs units. It is easy to verify that, for all t within the 
radiation era, these scales are comparable to the horizon 
size along the xy-plane (£h± = 2ci = 6 x 10 10 x t, in cgs 
units). In fact, 



0.67 £ H ± < Kj < 1-25 £ H ± 



(32) 



with the lower values (i.e. those for which X r j < £h±) 
corresponding to 7 < 0.5 and the higher ones to 0.5 < 
7 < 1. Confining ourselves to sub-horizonal scales (e.g. 
see [37]), we may conclude that, as regards favorism of 
gravito-magnetic instabilities in a cosmological model of 
the form (5), the values of 7 > 0.5 should be discarded. 

As regards the sub-horizonal perturbations, for a fixed 
k, Eqs (21) and/or (24) introduce to the evolution of 
the corresponding mode a critical time-scale {tk TJ ) in the 
following sense: As long as 



, ^ VA(7)/Q 3 (7) , 
1 < ( 7 + 1) fc - **- 



(33) 



the cosmological perturbation involved is unstable, i.e. 
n k < tk r j) < and grows exponentially with time [38] 

6p(t) , 8H(t) = A l exp (-1 J n k (t') d£ 



At 



x exp 



fP5~V P 5 ~ (7 + l) 2 Q 3 fc 2 i 2 



P 5 + VP* - (7 + i) 2 Q 3 fc 2 * 2 

(7)V^ - (7 + l) 2 Q 3 fc 2 * 2 



-ia( 7 )vTi 



where Ai stands either for A p or Ah, 

0(7) 



2( 7 + lV(22-9 7 )(3-7) 



(34) 



(35) 



and the exponentially decaying solutions - those with the 
plus sign in the exponents - arc neglected. With the aid 
of Eq (24) , Eq (34) is written in the form 



5p{k) , 6H(k) = At 



Wi-(^) 2 



l-(^7) 2 



^(7)^(7) 



(36) 



In the limit k —> k rj , Eq (36) yields 5p , SH = Ai (the 
constant amplitude of an oscillating perturbation), while, 
for k <^ k r j we obtain 



Sp(k) , 5H(k) > Ai ( ^ 



a(7)V-Ps(7) 



> A, 



(37) 



Accordingly, condensations of the corresponding length- 
scale (A rj ) will collapse in Jeans' fashion, while, at the 
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same time, the fluctuations of the magnetic field will 
grow, resulting in a gravito-magnetic instability, as well 
as in the amplification of the background magnetic field 
(in connection, see [25], [39]). 

The behavior of the (large scale) magnetized pertur- 
bations in the anisotropic model under consideration, is 
completely different than their behavior in the isotropic 
FRW model. In that case, both the Newtonian [11] and 
the GR study [13], [14] have verified a negative role of 
the magnetic field in the aggregation of matter on large 
scales. In particular, it has been shown that, on scales 
between the Jeans length and the horizon, the density 
perturbations may undergo a power-law evolution sim- 
ilar to that of the non-magnetized case (e.g. see [1]), 
but their growth-rate is reduced by an amount propor- 
tional to the magnetic field's strength. Therefore, in the 
isotropic case, the external field disfavors viable matter 
condensations. 

As regards the anisotropic case, the density perturba- 
tion at any spatial location can be obtained by a super- 
position of modes with different wave-numbers 



-0.4 



> 




FIG. 3: The spectral index v{pj) of the unstable perturbations 
as a function of 7, taking values in the range [| , 1). 



Sp(x, t) 



d 3 k 
(2tt) 3 



Sp(k) e 



(38) 



The strength of such a perturbation can be measured by 
the value of \Sp(x,t)\ 2 . From Eq (38), it follows that 
modes in the range (k , k + dk) contribute to \Sp(x,t)\ 2 
an amount proportional to \Sp(k)\ d k. Writing the con- 
tribution as 

\5p(k)\ 2 d 3 k = 4tt \6p(k)\ 2 k 2 dk = Air k 3 \5p(k)\ 2 d(\nk) 



we see that each logarithmic fc-interval contributes to 
\8p{x,t)\ 2 an amount V\ ~ k 3 \5p(k)\ 2 , which represents 
the so-called power spectrum [1]. With the aid of Eq 
(37) , the power spectrum of the MHD perturbations un- 
der consideration is written in the form 



r 3-2a( 7 )^P 5 ( 7 ) 



(39) 



which, for 7=4 reads Vk 



i reads Vk ~ k u y , while, in the limit- 
ing case 7 — » 1, it results in Vk ~ fc~ ' 6 . Notice that, 
although both cases correspond to the absence of the 
ambient magnetic field, the spectral indices involved are 
different. This difference arises due to the stiffness of the 
Universe matter-content as 7 — > 1 (c s — > c). The explicit 
form of the spectral index 



1/(7) 



k d 
Sp dk 



5p 



a(l)y/W) (40) 



as a function of 7 € [|, 1), is presented in Fig. 3. We note 
that, for every 7 in this range, we have f( 7 ) < 0. Once 
again, we verify that most of the power is concentrated 
on large scales (small wave- numbers) . 



IV. JEANS LENGTHS IN A RESISTIVE 
PLASMA 

It has been argued [2G] that, prior to recombination, 
the conductivity of the Universal plasma fluid is finite, 
although it may reach at very high values, of the order 
10 12 - 10 13 ml ohm). Accordingly, the study of MHD 
perturbations in a Bianchi Type I model has indicated 
that the finite conductivity (non-zero resistivity) could 
result in an exponential growth of the density perturba- 
tions. This argument was accordingly confirmed in the 
Newtonian limit [11], but not in the rclativistic one. In 
this Section, we shall attempt to re-examine the above 
statements in the rclativistic limit. The equations which 
govern the hydrodynamics of a magnetized, resistive fluid 
in the frame of GR, are identical to Eqs (1) and (2) being 
supplemented by the Maxwell equations 



2 

-OH* 



■u^+rjh^H^ (41) 



which coincide to Eqs (3), with the addition of an extra 
term due to the resistivity (77) of the cosmic medium. 

In this case, as regards the propagation of cosmological 
perturbations along the x-direction, Eq (16) reads 



(H 6H)[(2-~/)6+k 2 r)+in] 



1 + 7 



(l+ 7 )-m]<5p (42) 



Accordingly, the vanishing of the real part of the disper- 
sion relation results in 



22-97 
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(c 2 s + 
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2-7 
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Ri)9\, Su = 



(43) 
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Eq (43) can be cast in the form 
"22-9 7 



9 + kTr) 



k 777,4 + k 2 m,2 + mo 



(44) 



where, for every 7 £ [| , 1), 

J77.4 = C 2 S 7] > 

777 = {2- 1 )9(J 1 -R 3 )-9u 2 A R 1 
P 5 (7) 



18 (3-7) (l + 7) 3 i 3 



< 



and 



777 2 = C 2 S (2 - 7 ) 9 + 9 U\ + ( J X - R 3 ) 77 

0a(7) Ml) 



(45) 



(46) 



18(3-7)(l+7)r; 77 36(3 - 7 )(1 + 7 ) 2 



while 



fl 4 (7) = 2(l+7)(3- 7 )(27 7 2 + 367 + 16l) (47) 
+ (1 - 7) (37 - 1) (-157 2 + 174 7 - 931) > 

The allowed wave-numbers for stable (oscillating) per- 
turbations can be easily obtained by Eq (44), admitting 
that n 2 > 0. Accordingly, 



m4fc 4 + ra-ik 2 + 777o > 



(48) 



Since m\ — Am^rao > 0, the trinomial on the lhs of Eq 
(48) possesses the couple of real roots 



-m 2 1 



J (1 2) 27774 27774 



7772 — 477747770 



(49) 



For small values of 77 (77 — > 0), i.e. within the context 
of an iiieaZ-plasma cosmology, the corresponding wave- 
numbers behave as 
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S (7) 



± 
± 



36777 (3-7)(l + 7)t 

1 gaCr) 
36777 (3-7)(l + 7)t 

n(7)/Q 3 (7) 
(l + 7 ) 2 < 2 



(50) 



and therefore, now, the condition for unstable perturba- 
tions reads 



< 7i < k 2 < k 2 j2 



< k 2 < k 



A(7)/Q 3 (7) 
(1 + 7 ) 2 t 2 



(51) 



i.e. it is identical to the corresponding non-resistive one. 
On the other hand, in the limit of large 77, we obtain 



^(1.2) 72 7 (3-7)(l + 7)2*2 
£4(7) 



(52) 
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In accordance, as long as 
< k 2 < k 2 - 



^4(7) 



7 = 36 7 (3-7)(l + 7) 2 7;2 



(53) 



the condition (48) is not fulfilled and the corresponding 
cosmological fluctuations are unstable and grow exponen- 
tially with time (n 2 < 0). As a consequence, condensa- 
tions of the corresponding length-scale will collapse in the 
Jeans' fashion, revealing an indirect resistive instability 
[the rhs of Eq (53) is independent of 77] . 

In what follows, we compare the length-scales arising 
by the critical wave-numbers of the resistive and the non- 
resistive case. Our aim is to examine the contribution 
of resistivity in achieving more realistic (sub-horizonal) 
scales of gravito-magnetic instability, than those pre- 
dicted by the non-resistive fluid. To do so, we con- 
sider the ratio of the critical wave- numbers (53) and (24), 
which is written in the form 



k 2 

non 



£4(7) QM 
367(3-7)^5(7) 



(54) 



and is depicted in Fig. 4. First of all, we observe that 
for 7 = 5 and/or 7=1, i.e. in the absence of the am- 
bient magnetic field, we obtain k res — k non . This not an 
unexpected result, since 77 is introduced to the MHD anal- 
ysis through the Maxwell equation (41), which involves 
a non-vanishing magnetic-field strength. On the other 
hand, a large resistivity can modify quite significantly 
the scales involved in a gravito-magnetic instability and 
especially those around the value 7 = | (the peak in 
Fig. 4) , where the non-resistive analysis predicts a super- 
horizonal Jeans length [A non (7 = fi^) = 1-20 £u±(t)]. 
Accordingly, for 7 = |, we obtain 



Ar 



0.83 X r , 



< 



e H ±{t) (55) 



Therefore, we may conclude that, as compared to the 
ideal-plasma case, gravito-magnetic instabilities within a 
magnetized-fluid are particularly favored by a large re- 
sistivity, in the sense that it prolongs the range of the 
potential values of 7 (and therefore the number of the 
corresponding cosmological models) which admit sub- 
horizonal scales of such an instability (from ^ < 7 < | in 
the non-resistive case, to | < 7 < | in the fully-resistive 
one). 

This result, which is based on the fact that the critical 
(Jeans) wave-number in the resistive case is always larger 
than the corresponding non-resistive one (elf Fig. 4), has 
a clear physical interpretation: Indeed, the invariant form 
of Ohm's law for a locally neutral plasma is written as 



n 



(56) 



where, J M is the current density and is the Faraday 
tensor. According to it, besides the convective field, a 
non-zero resistivity favors also convective currents. As a 



8 



OJ 



eg 




0.98 



FIG. 4: The ratio kr e3 /k^ on of the unstable perturbations as 
a function of 7, taking values in the range [| , 1). 



consequence, for 77 + 1 0, the energy available to be ab- 
sorbed by the perturbation quantities is larger than that 
of the non-resistive case. This energy-excess will lead to 
the increase of the perturbations' frequency and hence to 
the amplification of their wave-number, as well (see also 
[25]). 



V. CONCLUSIONS 

Gravitational instabilities in a magnetized FRW Uni- 
verse, in which the magnetic field was assumed to be too 
weak to destroy the isotropy of the model, are known 
and have been studied in the past [13], [14], [16]. How- 
ever, mathematically speaking, the spatial isotropy of the 
FRW Universe is not compatible with the presence of 
large-scale magnetic fields. In fact, an anisotropic cosmo- 
logical model can and should be imposed for the treat- 
ment of magnetic fields whose coherent length is compa- 
rable to the horizon length [27]. Interesting homogeneous 
and anisotropic cosmological models of Bianchi Type I 
with an ambient magnetic field are known [28] and have 
been used to discuss astrophysical implications [6], [19], 
[25], [29], [31]. 

In this article, we have used the GR version of the 
MHD equations [18] in a homogeneous and anisotropi- 
cally expanding, self-gravitating, magnetized plasma, to 
study the gravito-magnetic excitation of cosmological 
perturbations. Accordingly, we have found that: 

• In a non-resistive medium: 

1. In contrast to the isotropic case, where the 
role of the magnetic field (as an agent oppos- 
ing the perturbations' growth) is clear, there 
is a wide class of anisotropic, magnetized cos- 
mological models - parameterized by 7 in the 



range ^ < 7 < 0.5, for which a gravito- 
magnetic instability can lead to sub-horizonal 
magnetized condensations [elf Eqs (29), (31) 
and (32)]. 

2. We have determined explicitly the power spec- 
trum of the unstable cosmological perturba- 
tions [Eq (37)], which indicates that most of 
the power is concentrated on large scales, close 
to the horizon [elf Eq(40)]. 

• In a resistive medium: 

1. The critical wave- numbers so obtained, ex- 
hibit a delicate dependence on resistivity [Eq 
(52)], resulting in the reduction of the corre- 
sponding Jeans lengths to smaller scales (well 
below the horizon) than the non-resistive ones. 
In fact, the non-zero resistivity induces con- 
vective currents, resulting in the enhancement 
of the perturbations' energy (frequency) and 
hence in the reduction of their comoving wave- 
length. 

2. Furthermore, the non-zero resistivity prolongs 
the range of the potential values of 7 (and 
therefore the number of the corresponding cos- 
mological models) which admit sub-horizonal 
scales of such an instability (from i < -v < 
in the non-resistive case, to | < 7 
fully-resistive one) [elf Eq (54)]. 



I < 7 < I 
3 — ' ~ 2 

< § in the 
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Appendix A 

The vanishing of the imaginary part, Di(k,n), of the 
dispersion relation (17) results in 
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1 + 7 
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where 



R2 = g (37 + 1) 



(A2) 
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and 



Accordingly, for n > we obtain k > kf , where 



Ra — — 



2(1 + 7) 
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(A3) 



In this case, the requirement n 2 > for oscillating 
(stable) fluctuations leads to the critical (Jeans) wave- 
number 



k 2 - 



1 + 7 



^ 9 
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kf 



1+7 



--(0-i? 4 )(2- 7 ) 



x{-J 1+ R 3 --(8-R 4 ) (2- 7 ), 



-Ji + R 3 



1 + 7 



[i?i + (1 + 7)^2] (A4) 



Once again, for every 7 <G [| , 1) the rhs of Eq (A4) 
is of the form (24), thus remaining positive for all t. 
Therefore, in a non-resistive fluid, the small-scale cos- 
mological perturbations with k < ki, arc unstable and 
grow exponentially with time. Hence, formations of the 
corresponding length-scale will collapse as a result of a 
gravito-magnetic instability. 

On the other hand, as regards the propagation of MHD 
fluctuations in a resistive fluid, the vanishing of the imag- 
inary part of the corresponding dispersion relation yields 



1 + 7 



[JZi + R 2 (1 + 7)] 



(A6) 



In this case, the critical wave-numbers so obtained 
are weakly-dependent on resistivity. 



(kf ), are weakly-dependent on resistivity. In particu- 
lar, around recombination the resistivity is related to the 
background temperature (T) according to [11], [40], [41] 



H 




= k 2 




l + 7_ 





1 + 7 9 



+J 1 -R 3 + -(e-R 4 )(2- 1 ) 6+ 



1 + 7 



[Rt +R 2 6(1 + 7)] 



(A5) 



V = Vn T- 3 ' 2 (AT) 
where r\ (in geometrical units) is measured in sec. At the 
late stages of the radiation epoch T = 2.7 (1 + z), where 
z is the cosmological redshift, while r/„ = 7 X 10 10 for a 
weakly ionized gas and r\ n = 10 13 for the fully ionized 
one. Therefore, at z rec ~ 1500, the product rj x 9 (in 
geometrical units) is dimensionless, admitting the value 



T) x 



2x 



10" 



(1+7) 



(for a weakly ionized gas). For this rea- 
son, the influence of resistivity in Eq (A6), as compared 
to Eq (A4), is of minor interest, since it simply adds a 
term of the order 10~ 8 to the evaluation of k 2 . 
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